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ABSIBAd 


The compurattonal aspects of modeling material failure In structural wood 
members are presented with particular reference to vector processing aspects. 
Wood members are considered to be highly orthotropic. Inhomogeneous, and 
discontinuous' due to the complex microstructure of wood material and the pres- 
ence of natural growth characteristics such as knots, cracks and cross grain 
In wood members. The simulation of strength behavior of wood members is 
accomplished through the use of a special purpose finite element/fracture 
mechanics routine, program STARW (Strength Analysis Routine for .Wood). Pro- 
gram STARW employs quadratic finite elements combined with singular crack tip 
elements In a finite element mesh which accounts for the complexities Inherent 
In wooa structural members. The need to use a highly refined finite element 
mesh to adequately model material behavior, results In the formulation of 
thousands of simultaneous equations which must be generated and solved repeat- 
edly to model the nonlinear failure process which occurs. The availability of 
the CYBER 205 at Colorado State University has made Implementation of program 
STARW at the level described not only possible, but also relatively economi- 
cal. Vector processing techniques are employed in mesh generation, stiffness 
matrix formation, simultaneous, equation solution, and material failure calcu- 
lations. The paper addresses these techniques along with the time and effort 
requirements needed to convert existing finite element code to a vectorized 
version. Comparisons In execution time between vectorized and nonvector I zed 
routines are provided. 
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INTRODUCTION 


Accurate knowledge ot the strength of a structural member Is essential 
Information to the design engineer concerned with structural safety and effi- 
cient material use. A means to predict material strength is necessary, since 
all materials exhibit some variability In strength and It Is not feasible to 
pnyslcal ly test every structural member to determine Its load carrying capa- 
city. The sophistication of strength prediction models have generally 
advanced, not only with the discovery and refinement of new computational 
methods, but also w h the Increase In computer capabilities which enable 
erflclent application of the new methods. 

In the case ot wood structural members, the current strength prediction 
method Is a highly approximate procedure based on empirical concepts from the 
1930's. This results In a strength prediction that Is relatively uncertain. 
The current strength prediction procedure Is based on the results of physical 
tests because until now It has not been possible to mathematically model wood 
member failure and rationally predict strength. The most obvious difficul- 
ties; orthotropic material properties, the presence of knots and associated 
grain deviations, and the presence of cracks from seasoning and partial 
material failure, can now be successfully modeled with program STARW (STrength 
Analysis fioutlne for Jtfood) (2). 

The nature of the nonlinear failure modeling process, presents a computa- 
tional problem of such a large magnitude that It can not be efficiently 
accomplished on computers that do not have the capacity of a CYBER 205. Pro- 
gram STARW represents a case where modest effort In Invoking vector processing 
syntax has not only made Implementation of the program possible, but has also 
resulted In a relatively economical solution. 
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AN OVERVIEW ON- MATHEMATICALLY MODELING WOOD ME MBER FAILU RE WITH STARW 


Program STARW uses two-dimensional ortholroplc finite elements to model 
behavior In the longitudinal -transverse plane of a loaded wood member. Ten- 
sile load Is applied In the longitudinal direction as shown In Fig. 1. 



Figure 1. Loaded Wood Structural Member 
(Long! tud Inal -Transverse Plane) 

A knot In a structural specimen of wood creates localized grain deviation 
as Indicated In Fig. 1. This grain deviation has an extremely Important 
effect on stress distributions at locations near the knot (3). An Iterative 
procedure to locate mesh coordinates corresponding to the grain deviation 
around a knot is employed In program STARW.- This procedure relates distortion 
of wood grain around a knot to streamlines of laminar fluid flow around an 
elliptical object and has therefore been named the "flow-grain analogy" (4). 
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Utilizing the flow-grain analogy, a representative finite element mesh Is 
automatically constructed of eight node quadrl lateral elements, six node tri- 
angular elements, and eight node singular elements. Since tangential elastic 
stiffness of wood may be as little as 1/20 of the longitudinal elastic stiff- 
ness, all three types of finite elements are required to model different elas- 
tic material behavior In the longitudinal and tangential directions. 
Appropriate elastic stiffness values for each element are automatically 
assigned. 

Singular elements are used to model material behavior around the tip of 
cracks that form as the load on the member Is Increased. These elements were 
developed using theory from linear elastic fracture mechanics (1). Experimen- 
tal Investigations have indicated that cracks In structural lumber will usu- 
ally form ana propagate along a grain line.- Thus, cracks are modeled by pro- 
gram STARW by "unzipping" the finite element mesh along the material separa- 
tion ana placing the singular elements arouna the crack tip. A resulting fin- 
ite element mesh Is shown In Fig. 2. The "unzipping" process and placement of 
tne singular elements are pertormed automatically upon cue by the user when 
the appropriate failure conditions are Indicated In the program output. 
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Direction of Applied Stress 




Around the Crock Tip 


Figure 2. Example- Finite Element Mesh Including Crack 


The output directly calculated from each analysis Is as follows: 

1) Horizontal and vertical displacement at each node In the mesh. 

2) Stresses for each element, paral I e'l-to-grai n, perpendlcu lai — to- 
graln, and shear. 

3) Stress Intensity factors resulting from the use of singular ele- 
ments. 


4) A failure summary that Indicates to the user what appropriate 
action should be taken to model the next step In the failure pro- 
cess. 
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The stress Intensity factors directly reflect the strength of the stress field 
around the crack tip. The stress Intensity factors are compared within the 
program to a fracture criteria for structural wood members to determine If the 
existing crack propagates at a given applied load. The element stresses are 
compared to a failure criteria for structural wood members to determine If a 
crack will form near the element under consideration. The results of these 
comparisons are expressed In the failure summary. 

Analyses are performed repeatedly with stress and stress Intensity fac- 
tors monitored at each step and compared within the program logic to the 
fracture/fal I ure criteria. As the load on the member Is increased, more 
cracking and material failure occurs. The user, based on the information In 
the failure summary and the overall stress picture, gives the program the 
necessary Information to model the successive step In the failure process. In 
the future, as research progresses, program logic will be expanded to Include 
the decision making process the user currently makes based on the failure sum- 
mary. Failure may be continually modeled In this fashion until the member 
unuer consideration has failed to the point where It cannot resist an Increase 
In load. At this point, the predicted strength Is realized. In studying the 
behavior of a wood member, 30 analyses may typically be performed before the 
member reaches Its capacity. A simplified diagram of the failure model Is 
contained In Fig. 3. 
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Figure 3. Strength Prediction Model 


ASPECTS AND IMPLICATIONS OF VECTQRLZATIflN 

For each analysis, program STARW performs five general sets of computa- 
tions: 

1) Generation of a suitable finite element mesh using the flow-grain anal- 
ogy and an unzipping process to Include cracks. 

2) Formation of a set of simultaneous equations which may be 2000 to 5000 
equations In length. 

3) Solution of the simultaneous equations using Gauss elimination. 
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4) Calculation and coordinate transformation of element stresses based on 
the solution vector and the element grain angles. 

5) Computations with the fal lure/ fracture criteria using element stresses 
and stress Intensity factors as Input. 

Routines Included In Items 1 through 4 existed In limited form and were 
executed for small problems on a CYBER 720 prior to application on the CYBER 
205. Failure calculations In Item 5 and additional mesh generation capabili- 
ties were added and designed specifically for use on the CYBER 205. After 
compiler Induced vectorlzatlon proved to be Inadequate, In significantly 
reducing execution time. It became apparent that it was essential to expli- 
citly vectorize selected portions of the existing routines. At the same time. 
It was not the primary goal of the project to expend unlimited effort to 
achieve the maximum In vector processing, rather the goal was to produce a 
powerful research tool that could be economically Implemented. The bulk of 
the conversion (and execution time savings) were achieved with modest effort 
after becoming familiar with vector processing syntax. 

To date, a means to vectorize the Iterative solution of the fluid mechan- 
ics equations contained In the flow-grain analogy has not been established. 
This Is not of great concern since, as in many finite element routines, mesh 
generation does not account for a significant portion of the total execution 
time. However, the unzipping of the finite element mesh to model cracks 
Involves, In part, a uniform renumbering of nodal points. This renumbering Is 
easily accomplished with basic vector commands since nodal coordinates are 
stored In vector form. 
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Formation ot the set of simultaneous equations can typically take from 5 
to 50 per cent of total execution time In a unvectorized finite element 
analysis. In program STARW, a 16 by 16 element stiffness matrix must be con- 
structed for each element and properly combined with other element stiffness 
matrices to form the coefficient matrix (global stiffness matrix) of the 
simultaneous equations. Formation of the 16 by 16 matrix Involves dot pro- 
ducts or vectors of length 16. Some time savings Is attained here through the 
use of the CYBER Q8SD0T command even though the vector length Is rather small. 

Solution of the simultaneous equations typically requires 40 to 90 per- 
cent of the total execution time of a finite element analysis. The 90 percent 
figure Is not uncommon for large two-dimensional analyses. Therefore, large 
time savings can be attained by vectorizing the solution algorithm alone. In 
program STARW, Gauss elimination is used to. decompose the global stiffness 
matrix, followed by a back substitution to obtain the solution. For the prob- 
lem under consideration the stiffness matrix Is banded and symmetric, and 
therefore, only the upper diagonal half of the matrix Is stored. Furthermore, 
If the global stiffness matrix Is stored In columns rather than rows, then 
adjacent terms In a row of the global stiffness matrix will be stored contigu- 
ously. Since Gauss elimination Involves operations of one row upon another, 
by storing the matrix as described, each row will be a vector. "Gather" and 
"scatter" vector formation commands are unnecessary. Gauss elimination 
Involves operations on the matrix rows In a number of nested DO loops. Vector- 
Ization of even the Inner most loop results In large time savings. Back sub- 
stitution Involves repeated dot products of previously formed vectors. This 
can again be easily accomp 1 1 shed w I th the CYBER Q8SD0T command. An unvector- 
Ized and otherwise Identical vectorized portion of the back substitution Is 
shown In Fig. 4 to Illustrate typical vector I z at I on. 
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DO 460 J * 2, JEND 

J1 « 11 + J - 1 

B( 1 1 ) = B( i 1 ) - A(J,m * B(J1) 

460 CONTINUE 

LE * JEND - 1 

J1 = 11 + 1 

B( 11) * B(I1) - Q8SD0T <A<2, II; LE), B(J1; LE)) 

Figure 4. Example DO Loop and Corresponding Vector Syntax 

With the solution of the equations established, element strains and 
stresses can be calculated In global coordinates. Since this calculation Is 
essentially the same for every element, and care Is taken to store the neces- 
sary quantities In vector form, basic vector operations accomplish this task. 
The solution vector Is found In the global coordinate system and thus the cal- 
culated stresses are also expressed In this system. It Is desfreable, however, 
to know the stresses In the coordinate system of each element or the 
perpendlcular-to-graln and paral lel-to-graln directions. The element stresses 
must be transformed according to the element grain angle. Since the element 
grain angles are stored contiguously and In order, this computation can be 
accomplished with basic vector commands. 

To complete an analysis, the stresses and stress Intensity factors for 
cracks must be Inserted Into the fal I ure/ fracture criteria. The 
fal lure/fracture criteria Interfaces the mathematical results from an analysis 
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+o the real life failure actions. Required Information Includes the maximum 
stresses and their locations within the flow-grain mesh. Since stresses are 
stored In element order In vectors, this Information can be obtained much 
quicker and more easily by using CYBER Q8 commands than with scalar search 
al gorl thms. 

To put the vector I zat I on discussed Into perspective, a typical problem 
was analyzed using unvectorized and vectorized routines. Since unvectorlzed 
versions of the mesh generator (Item #1) and the maximum stress searching rou- 
tine (Item #5) do not exist, vectorized routines had to be used for both sides 
of the example. The example problem consisted of 4180 degrees of freedom 
(equations) and for simplification no cracks were Included. The corresponding 
CPU execution times for different phases of the analysis are shown In Table 1. 


Table 1. Efficiency of Execution Time for Vectorized Routines 



Unvectorized 
time in sec. 

Vectorized 
time in sec. 

Efficiency 

unvect/vect 

Mesh Generation 

1.90 

1.90 

1.00 

Stiffness Matrix Formation 

4.84 

2.80 

1.73 

Solution of Equations 

97.87 

4.91 

19.90 

Ml SCELLANEOUS COMPUTAT 1 ONS 

5.05 

4.60 

1.10 

Total 

109.66 

14.21 

7.70 


As clearly shown for this problem, the vectorized equation solver was 20 
times faster than Its otherwise Identical unvectorlzed version. This savings, 
along with other vector I zat l on, reduced analysis time by nearly a factor of 
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eight. One will note that while the miscellaneous computations were somewhat 
Insignificant In the unvectorized analysis, they take on new importance In 
tne vectorized analysis. Additional effort may be well spent in further vec- 
torlzatlon of the miscellaneous computations. 

CONCLUSIONS 

Failure In wood members Is being successfully modeled and analytically 
Investigated In greater detail than before possible through Implementation of 
program STARW on the CYBER 205 (2). An understanding of material failure Is 
essential to accurately predict member strength and to safely and efficiently 
use tne material in engineering application. 

Vectorlzatlon of program STARW has reduced an unwleldly and expensive, 
nonlinear failure modeling method Into an efficient research tool. Vectorlza- 
tlon of existing routines need not be a lengthy and laborious effort to 
achieve execution time savings. It has been shown that careful organization of 
operands Into vectors and modest effort In invoking vector syntax can cut pro- 
gram execution time by a factor of nearly 8 for a typical problem In this 
research. The largest savings Is realized In the solution of the simultaneous 
equations. 

While use of program STARW Is expected to provide new Information on 
fracture and failure In wood members, the availability of machines with the 
capabi lines of the CYBER 205, in general holds promise for advances In the 
analytical modeling of all materials. These advances In research will ini- 
tiate new applications of materials and more efficient and reliable use of 
materials In existing applications. 
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